Resemblances and differences in mechanisms of noise-induced resonance 
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Systems showing stochastic resonance (SR) or coherent resonance (CR) share some features, in 
particular the nearby periodic character of the signal. We show that in spite of this resemblance 
the different underlying dynamics can be detected in experimental data by studying the histogram 
of inter-spikes times and some statistical properties like two-times correlation functions. We discuss 
the possible relevance for climate modeling. 
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I. INTRODUCTION 

The mechanism of stochastic resonance (SR) was ini- 
tially introduced as a possible explanation of climate 
changes on long time-scales pj. During the last two 
decades it has been applied to a wide class of systems 
ranging from analog circuits, neurobiology, ring lasers, 
systems with colored noise, etc; for a review see Ref. 0- 

The prototypical system showing SR, which is also 
the original one used to model climate changes, is the 
stochastic differential equation 



dx 
~dl 



dV(x,t) 
dx 



(1) 



where 77 is a Gaussian, white noise with (77(f)) = and 
(ji{t)ri(t')) = S(t—t'), D measures the noise intensity and 
V(x, f) a is double well potential with a time periodic 
term 

4 2 

V(x,t) = — - — + Arcos(27rf/T) . (2) 

In the case of a stationary potential, i.e., ^4 = 0, the 
jumps between the two minima at x = — 1 and x — 1 
are independent events whose probability distribution is 
approximately Poissonian Q. Using simple arguments 
based on the Kramers exit-time formula |4|, it can be 
shown that there is range of values of D, T and A where 
SR is present, i.e., the jumps between the two minima 
(close to —1 and +1 if A is sufficiently small) are strongly 
synchronized with the forcing and that the probability 
distribution function (PDF) of the jumping time r has a 
relatively sharp peak around T [J, • 

The phenomenon of SR provides one example of the 
nontrivial role that noise can play in dynamical systems 
with an external periodic forcing. Besides SR, there exist 



other examples of the "constructive role" of noise, e.g., 
one can have a synchronization of trajectories generated 
by different initial conditions and the same noise real- 
ization [5j. Our interest will focus on cases where noise 
can enhance periodic behavior, e.g., the so-called coher- 
ent resonance (CR) and the noise-induced dynamics in 
systems with time delay (ND). 

The phenomenon of CR [|| has been found in mod- 
els describing excitable systems that occur in different 
fields like chemical reactions, neuronal and other biologi- 
cal processes 0, 13 • The prototypical stochastic differen- 
tial equation used in this case is the FitzHugh-Nagumo 
system defined by: 



dx 

~db 

dy 

dt 



x" 
x + a + 



(3) 
(4) 



with e <C 1 so that the time evolution of x is much faster 
than that of y. For \a\ > 1 there is a stable fixed point, 
for \a\ < 1 there is an unstable fixed point and a limit 
cycle. The cycle consists of two pieces of slow motion 
connected by a fast jump. If \a\ is slightly larger than 1 
the system is excitable ((J , i.e., small deviations from 
the fixed point may generate large pulses (also called 
spikes) 0- Moreover, in this case, one finds that there is 
a range of values of the noise intensity D such that CR 
appears, i.e., roughly periodic noise-excited oscillations 
are present, resembling the SR oscillations [||. 

The prototypical example for ND 01 is the over- 
damped particle motion in the double-well potential 
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FIG. 1: PDF of the inter-spikes time in the system described 
by Eq. Q showing stochastic resonance. The external peri- 
odic force has period T = 100 and amplitude A — —0.15. 
The noise intensity is at the corresponding optimal value 
D* = 0.10. 



FIG. 2: PDF of the inter-spikes time in the system described 
by Eq. ijjy. The amplitude of the delay term is A = —0.15, 
and the delay time is T = 100. The noise intensity is at the 
corresponding optimal value D* = 0.10. 



V(x(t),x(t-T)): 

dx(t) dV{x{t),x{t-T)) 



dt dx(t) 

= x(t) - x(t) 3 - Ax(t — T) + V2D-q (5) 

where T is the delay. It is not very difficult to realize that 
the delay term Ax(t — T) has a role similar to that of the 
periodic forcing in Eq. , |2J . Accordingly, in a certain 
range of parameters' values , there is a sort of periodic 
motion with period T or 2T (this depends upon the sign 
of A). This ND equation had been proposed as a model 
for some climate changes |Tlj . 

Although SR, CR and ND, are similar phenomena, in 
the sense that their time evolution is nearly periodic, 
there are also some important differences. For example: 
a) due to the presence of the term Ax(t — T), Eq. is 
in fact an infinite dimensional system since in order to 
determine x(t) for t > one has to specify x(t') with 
— T <t'<0. On the contrary for Eq. it is sufficient 
to know x(0). b) in the case of SR the periodicity is due 
to the external forcing while in the CR case the period- 
icity has an internal origin, i.e., the periodic motion is 
due to the intrinsic dynamics and, at variance with SR, 
its period cannot be changed by tuning external control 
parameters. This difference can play an important role 
in, e.g., the context of climate changes and glaciation, for 
more details, refer to the last Section. 

The aim of this paper is to analyse the differences 
among SR, ND and CR and their possible relevance to 
applications. In Sect II we will briefly review some prop- 
erties of the PDF of the inter-spikes times for SR, ND and 



CR. In particular, we recover a recent result [l3| showing 
that CR and SR are not conflicting or excluding mecha- 
nisms, i.e., the same periodically driven system, e.g., the 
one given by Eqs. Il.''>l 111 with a time-dependent param- 
eter a(t) = ciq + aicos(2irt/T) can present a transition 
from SR to CR behavior when the noise intensity D is 
increased . In Sect III we show that, in spite of some 
resemblance, SR, ND and CR exhibit different statisti- 
cal features which, at least in principle, can be detected 
in experimental data. In particular we have that for SR 
the correlation function C(r), after a transient period, is 
periodic and does not relax to zero. On the contrary, for 
the CR and ND cases C(r) shows damped oscillations. 
Sect IV is devoted to general remarks and conclusions. 
In particular we deal with the potential relevance of the 
differences between SR, CR and ND to climate modeling. 



II. STATISTICS OF INTER-SPIKES TIMES 

One basic feature shared by the three models intro- 
duced in the previous Section, is the presence of two char- 
acteristic states: two equilibria in the case of SR and ND, 
a rest state and an excited one in the FitzHugh-Nagumo 
model. Jumps between these states are made possible 
by the noise. Moreover, it turns out that there exists an 
optimal value of the noise intensity such that this jump- 
ing becomes approximately periodic, i.e., the typical time 
between two consecutive transitions is roughly constant. 

According to the terminology of biological systems, 
where CR was originally introduced, we will refer to 
the time interval between consecutive transitions as the 
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FIG. 3: PDF of the inter-spikes time in the FitzHugh- 
Nagumo model Eqs. (l.'il lfi with parameters' values ao = 1.05, 
ai = 0.04, e = 0.01 and T = 10 > T* = 3.7 for differ- 
ent noise intensities. From top to bottom, D — 2.5 x 10~ 5 , 
D = 4.3 x 10 -5 , D = 1.1 x 10~ 4 , D = 2.3 x 10~ 4 and 
D = 4.3 x 10~ s . 

"inter-spikes time r". In the case of the system (!.''>! II) 
where the dynamical variables show well defined maxima, 
the definition of the inter-spikes time r is rather natural. 
In the case of the other two systems 111 121) and © one can 
define r as r = t n+ \ — t n where t n is the n— th crossing 
time, i.e., x(t n ) = and x(t n ) > 0. A measure of the 
signal's periodic character is provided by the normalized 
variance NV of inter-spikes times, NV = y/Var(r) / (r) . 
For generic noise intensity the transitions occur at ran- 
dom times, and p(r), the PDF of inter-spikes times, is 
weakly localized, i.e., NV is of order 1. We define the 
resonant or optimal value of the noise intensity, and de- 
note it by D*, as the value of D for which the system has 
minimal normalized variance NV. 

In Figure ^ we show the PDF of the inter-spikes time 
interval in the case of SR for the system described by 
Eq. at the optimal noise intensity D* = 0.1. We 

see that p(r) is peaked around r = nT with n — 1, 2, , 

and that the envelope is approximately exponential. This 
feature can be easily explained as follows. Consider the 
trajectories x(i), starting at t — from the favored well, 
i.e. x(Q) = 1 if A < 0. In the case of SR, at t close to T/4 
many of the trajectories will jump onto the other mini- 
mum at x = — 1 and after half period they will jump back 
again onto x = 1. However, a fraction of the trajectories 
remains in the "wrong position" (i.e. in the unfavored 
well) for t close to T. Calling P the probability of this 
event, we have that the integral of p(r) around T, say for 
r G [0.5T, 1.5T], is (1 - P). The events with inter-spikes 
time r ~ nT correspond to trajectories x(t) which are in 



FIG. 4: As in Fig. (0) but with T = 1< T* = 3.7. 



the the "wrong" minimum at T, 2T, nT. Taking into 
account the periodicity of the forcing and assuming that 
the system's memory is much shorter than the external 
period T we have that the probability to have x(t) in the 
"wrong" minimum at t ~ kT under the condition that 
x(t) was in the the "wrong" minimum at t ~ (k — 1)T, 
does not depend on the behavior for t < (k — l)T. There- 
fore the integral of p(r) around nT is (1 — P)P n ~~ 1 ~ e~ cn 
with c = — In P. The envelop of the inter-spikes-interval 
histogram has recently been computed by Berglund and 
Gentz [l3| in a more general and rigorous setting. 

In the case of CR (not shown) and ND, see Figure |2 
the PDF of the inter-spikes time interval is peaked around 
T where T is the characteristic internal time of the sys- 
tem, namely the delay time in the ND case, and the pe- 
riod of the limit cycle in the CR case. The parameters of 
the memory term in the ND-model (0 have been chosen 
in order to emphasize the similarities with the SR-model, 
i.e., the amplitude A and delay time T are identical to the 
amplitude and period of the external forcing in Eq. lfT|l. 
One consequence of this choice is that the intensity of 
the optimal noise coincides with that in the SR-modcl. 
Notice that, at variance with the SR case, no peaks are 
present at multiples of T. 

In summary, in all three systems there exists an op- 
timal noise intensity which produces a roughly periodic 
signal x(t). At this intensity a sharp peak appears in 
the inter-spikes time PDF. In the case of CR and ND 
the PDF has only one maximum at r = T, On the other 
hand, in the case of SR other maxima appear at nT. This 
effect can be considered as one of the distinctive marks 
of the stochastic resonance. In contraposition to the CR 
and ND cases, the periodicity of the signal in the SR case 
is induced by an external periodic force, which triggers 
the jumps at fixed times t = nT in such a way that the 
system synchronizes with this external "clock" . 

The combined effects of noise and nonlinearity can re- 
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FIG. 5: Conditional average (a;(t)|a;o) (solid line) and condi- 
tional variance a 2 (t) (dashed line) in the case of SR for the 
system described by Eq. The initial position is xq — 1, 
i.e., in the favored well. The external force has period T = 100 
and amplitude A — —0.15. The noise intensity is at its op- 
timal level D — 0.10. The average is taken over N = 10 4 
realizations of the noise. 



suit in more complex behavior |lj, 1 1 5| . In particular it 
has recently been shown in [l^ that SR and CR can co- 
exist in the same system. This behavior occurs, e.g., in 
the model II) when small oscillations are imposed on 
the control parameter a(t) = ao + ai cos(2nt/T). We con- 
sider only clq and oscillation amplitudes a\ such that the 
control parameter a(t) never crosses the critical value, 
i.e. a(t) > l,Vt. The presence of these small oscilla- 
tions determines privileged times t n — (2n + l)T/2 at 
which the system is closer to the excited state, i.e., a(t n ) 
approaches 1 from above, and a noise induced transi- 
tion is facilitated. Analogously to the case of the double 
well system with periodic forcing, there exists an optimal 
noise intensity for which a regular, quasi-periodic behav- 
ior emerges. More precisely: at low noise intensity, p(r) 
has rather sharp peaks at r = nT with an approximately 
exponential envelope, showing all the features of SR, see 
the two upper panels in Figure At higher noise inten- 
sity, as in the two lower panels in Figure [31 the small os- 
cillations in a(t) become irrelevant and the system shows 
CR as if the control parameter were fixed at its mean 
value arj. The shape of p(r) behaves accordingly as the 
noise increases: the maxima of SR diminish and a single 
peak with an exponential tail appears at r = 7 1 * where 
T* is the period of the system's limit cycle. 

Moreover, when the external force period T is shorter 
than the internal, limit-cycle period T* there are peaks 
only for nT > T* because once the excited state is 
reached the system needs at least a time T* in order to 
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FIG. 6: Same as in Fig. I^Jin the case of CR for the system 
described by Eqs. 1MIH with a — 1.05, e = 0.01 and optimal 
noise intensity D = 2.5 x 10~ 3 . The initial position is xo = 
2, yo = 0.8, i.e., in the excited state. The average is taken 
over N = 10 4 realizations of the noise. 



relax back into the excitable, rest state and restart the 
whole cycle (see fig-EJ. 

III. CONDITIONAL AVERAGES AND 
CORRELATION FUNCTIONS 

If one would observe just a single trajectory, SR, CR 
and ND would appear rather similar since the three cases 
would present us with a nearly periodic x(t). An analysis 
based on Fourier spectra, as it is often done, would rein- 
force this picture. In the previous Section we contrasted 
the multi-peaked PDF of inter-spikes time in the SR case 
with the one-peaked PDF in the CR and ND cases. In 
this Section we bring to the fore some statistical proper- 
ties which are present in the SR case but are absent both 
in the CR and in the ND cases. 

Consider an ensemble of N trajectories {x^(t),n — 
l,....,N 1}, sharing the same initial conditions 
a;( n )(0) = xq, but with different realizations of the noise 
r](t) and compute from Eqs. (JIJ, (l.'il ll) and (J5J the condi- 
tional average (x(t)\xo), 

1 N 

(x(t)\x ) = -J2 x{n) (t), (6) 

n=l 

and the conditional variance 

<7 2 (t) = (X 2 (t)\x )-(x(t)\x ) 2 
1 N 

= ^& (n) (*)-(z(*)M] 2 , (7) 

71=1 



5 



o 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 




1 2 3 4 5 6 7 
t/T 



9 10 



O 




FIG. 7: The periodic correlation function C(r) of the SR 
model Eq. Parameters' values as in Fig.Q 



In Figs. Eland El we show (x(t)\x ) and cr 2 (t) as func- 
tions of the time t, for the SR case and the CR case re- 
spectively. In both cases the noise-intensity values D are 
the optimal ones. In the SR case (x(t)\xo) does not relax 
to zero at large times t and, ignoring the initial transient, 
it is periodic. The conditional variance <J 2 (i) reaches its 
minima when the absolute value of the conditional av- 
erage (x(t)\xo) reaches its maxima. The largest values 
for <J 2 (t) are achieved when (x(t)\xo) is around zero. In 
other words, the main uncertainty in the process occurs 
around T/4, 3T/4, 5T/4 and so on, i.e., when the jumps 
between the two minima take place. 

As it can be seen in Fig El the behavior found in the 
CR case is different. Even with the noise intensity at its 
optimal value, after a few damped oscillations both the 
conditional average and variance relax to the constant 
values ((x(t))) and ((x 2 {t))) - ((x(t))) 2 , where ((•••)) 
indicates a time average. 

This behavior underlines the intrinsic difference be- 
tween the SR and CR mechanisms. In the case of SR, 
the presence of an external synchronizing force, makes 
the transitions from state +1 to state —1 to occur around 
preferred times. A set of independent replicas, initially 
localized in one of the two wells, will therefore quickly 
reach a periodic configuration, with the maximum prob- 
ability localized in the time-dependent favored well. On 
the contrary, in the CR case, there are no externally de- 
fined preferred times for the transitions, therefore each 
replica quickly loses its initial synchronization with the 
other ones and after a few periods the jumps occur at 
different times for different replicas. 

The correlation function 
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FIG. 8: Correlation function C(t) of the y variable for the 
FitzHugh-Nagumo model described by Eqs. (0 , EJ • Parame- 
ters' values as in Fig. El 



to zero while in the SR case it remains periodic with non 
decreasing amplitude, see Figs. (0 and JHJ. If one defines 
a correlation time t, as 



C(r) 2 dT, 



(9) 



behaves similarly to (x(t)\xo): in the CR case it relaxes 



one finds that r c diverges in the SR case while it remains 
finite in the CR case and that, as a function of the noise 
intensity, it attains its maximum r c = 1.03 at the optimal 
noise intensity value. The behavior in the ND case (not 
shown) is qualitatively very similar to that obtained in 
the CR case. 



IV. CONCLUSIONS AND DISCUSSION 

The numerical study of SR, CR and ND presented in 
this paper is focused on contrasting statistical features 
which would be difficult or impossible to detect by sim- 
ply looking at the spectrum since the spectrum shows a 
peak that achieves a maximal sharpness for the optimally 
chosen noise intensities in all the cases we have studied. 
As we have seen, it is possible to distinguish between the 
different underlaying dynamical mechanisms by studying 
the PDF of inter-spikes times, conditional averages and 
time-delayed correlation functions. 

We believe that these results are of interest not only in 
the context of dynamical systems but also in the study 
of certain climate phenomena. For time scales of order 
O(10 5 ys) and larger, Milankovich has proposed that 
Earth's climate has been determined by the influx of so- 
lar energy to such an extent that the fluctuations in, e.g., 
the global mean temperature and seasonality must have 
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been closely correlated with the variations in the incom- 
ing energy flux due to the periodicities in Earth's orbit. 

Another possible issue, for which the presented results 
are potentially interesting, is the so-called Dansgaard- 
Oeschger(DO) events [17J which have been inferred from 
the study of Late-Pleistocene ice cores and marine sedi- 
ments. These measurements show rapid warmings of the 
atmosphere followed by a much slower decay back into 
the average glacial conditions. The warmings took place 
on a time scale of a few decades while the relaxation 
back into glacial temperatures lasted some centuries up 
to millennia. They seem to have occurred at intervals of 
1.500 ± 200 years or integer multiples hereof They 
were absent during the Holocene, i.e., during the 10 4 
years before present. The discovery of these rapid warm- 
ings led to proposing a number of possible explanations, 
some of them favoring the internal origin of the period 
approximately equal to 1.500 years 0, |2(], while 
other explanations assume a similar period due to an ex- 
ternal astronomical forcing [22^. i.e. a SR-type scenario. 
In particular, it was shown that in an ocean-circulation 
box model and within an appropriate parameters' range, 
the purely deterministic system has a fixed point and 
does not show any time dependence while the addition 
of noise leads to the generation of spikes with a well de- 
fined inter-spikes time interval |23| . i.e., that coherence 
resonance is present in this ocean circulation model. As 



discussed in Section II, now we know that, at least in 
some systems, it is possible to observe either CR or SR 
behavior depending upon the noise intensity. 

We have shown that it is possible to distinguish be- 
tween the SR and CR, e.g., by looking at time-delayed 
correlation functions and at the PDF of the inter-spikes 
times. Needless to say, in order to compute such cor- 
relation functions, or the PDF, a sufficiently long and 
accurate series of measurements is required. From the 
limited information about the DO events that has been 
extracted from the geological record it is difficult to de- 
cide in favor of one scenario or the other. Indeed, the 
PDF for the inter-spikes time is qualitatively in agree- 
ment with the one observed in the SR case. On the other 
hand the time values at which one observes peaks of the 
PDF do not correspond to known astronomical periods, 
accordingly, a CR scenario would seem more appropriate. 
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